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Abstract. Although various cosmological observations congruously suggest that our universe is dominated by two dark com- 
ponents, the cold dark matter without pressure and the dark energy with negative pressure, the nature and origin of these 
components is yet unknow. The generalized Chaplygin gas (gCg), parametrized by an equation of state, p = -A/p'Jj,^, was 
recently proposed to be a candidate of the unified dark matter/energy (UDME) scenarios. In this work, we investigate some 
observational constraints on it. We mainly focus our attention on the constraints from recent measurements of the X-ray gas 
mass fractions in clusters of galaxies published by Allen et al. (2002,2003) and the dimensionless coordinate distances to type 
la supernovae and Fanaroff-Riley type lib radio galaxies compiled by Daly and Djorgovski (2003). We obtain the confidence 
region on the two parameters fully characterizing gCg, A , = A/p^^^^g and a, from a combined analysis of these databases, where 
PgCgo is the energy density of gCg at present. It is found that A, = O.TO^^jJ jy and a = -0.09^q33, at a 95% confidence level, 
which is consistent within the errors with the standard dark matter + dark energy model, i.e., the case of o- = 0. Particularly, the 
standard Chaplygin gas (a = I) is ruled out as a feasible UDME by the data at a 99% confidence level. 

Key words, cosmological parameters — cosmology: theory — distance scale — supernovae: general — radio galaxies: general 
— X-ray: galaxies:clusters. 



1. Introduction 

Two dark components are invoked to explain the current cos- 
mological measurements: the cold dark matter (CDM) with- 
out pressure and the dark energy (DE) with negative pressure 
(for a recent review, see Peebles & Ratra 2003). The first one 
contributes f2,„ ~ 0.3, and is mainly motivated to interprete 
galactic rotation curves and large scale structure formation 
(e.g. Longair 1998), while the second one (Qde ~ 0.7) pro- 
vide a mechanism for acceleration discovered by distant type 
la supernovae (SNela) observations (Perlmutter et al. 1998, 
1999; Riess et al. 1998, 2001), and offset the deficiency of a 
flat universe, favoured by the measurements of the anisotropy 
of CMB (de Bernardis et al. 2000; Balbi et al. 2000, Durrer 
et al. 2003; Bennett et al. 2003; Melchiorri & Odman 2003; 
Spergel et al. 2003), but with a subcritical matter density pa- 
rameter f2,„ ~ 0.3, obtained from dynamical estimates or X- 
ray and lensing observations of clusters of galaxies(for a re- 
cent summary, see Turner 2002). There are a huge number 
of candidates for DE in the literature, such as a cosmologi- 
cal constant A (Carroll et al. 1992; Krauss and Turner 1995; 
Zhu 1998; Sahni 2002; Padmanabhan 2003), the so-cafled "X- 



matter" (Turner and White 1997; Zhu, Fujimoto and Tatsumi 
2001; Lima and Alcaniz 2002; Lima, Cunha and Alcaniz 2003; 
Gong 2004; Chen 2004), and quintessence (Ratra and Peebles 
1988; Caldwell et al. 1998; Sahni and Wang 2000; Gong 2002; 
Sahni et al. 2003; Padmanabhan and Choudhury 2003) etc.. 
However, neither CDM nor DE has laboratory evidence for 
its existence directly. In this sense, our cosmology depends on 
two untested entities. It would be nice if a unified dark mat- 
ter/energy (UDME) scenario can be found in which these two 
dark components are different manifestations of a single fluid 
(Padmanabhan & Choudhury 2002; Wetterich 2002; Matos & 
Urena-Lopez 2000). 

Recently, the generaUzed Chaplygin gas (gCg) was pro- 
posed as such a unification, which is an exotic fluid with the 
equation of state as follows 



PgCg = -A/p"cg 



(1) 



where A and a are two parameters to be determined. It was 
originally suggested by Kamenshchik et al. (2001) with a - 1, 
and later on extended by Bento et al. (2002) to gCg. This sim- 
ple and elegant model smoothly interpolates between a non- 
relativistic matter nhase (d = 0) and a neeative-oressure dark 
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2002) and moreover it admits a well established brane inter- 
pretation (Kamenshchik et al. 2001; Bento et al. 2002; Bilic et 
al 2002). It is promising to be a candidate of the UDME. Such 
a possibility has triggered, quite recently, a wave of interest 
aiming to constrain the gCg model using various cosmological 
observations, such as SNela (Fabris et al. 2002; Makler et al. 
2003a,b; Bean and Dore 2003; CoUstete et al. 2003; Silva and 
Bertolami 2003; Cunha, Alcaniz and Lima 2004; Bertolami 
et al. 2004), the CMB anisotropy measurements (Bento et al. 
2003a,b; Bean and Dore 2003 ; Amendola et al. 2003), the grav- 
itational lensing surveys (Dev, Jain and Alcaniz 2003, 2004; 
Silva and Bertolami 2003; Makler et al. 2003b; Chen 2003a,b), 
the X-ray gas mass fraction of clusters (Cunha, Alcaniz and 
Lima 2004; Makler et al. 2003b), the large scale structure (Bilic 
et al. 2002; Bean and Dore 2003; Multamaki et al. 2004), and 
the age measurements of high-z objects (Alcaniz, Jain and Dev 

2003) . But the results are disperse and somewhat controversial, 
with some of them claiming good agreement between data and 
the gCg model while the rest ruUng it out as a feasible UDME. 

In this work, we shall consider the observational constraints 
on the parameter space of gCg arising from the X-ray gas mass 
fractions of clusters of galaxies pubUshed by Allen et al. (2002, 
2003) and the dimensionless coordinate distances to SNela and 
Fanaroff-Riley type lib (FRIIb) radio galaxies compiled by 
Daly and Djorgovski (2003). We perform a combined analy- 
sis of these databases and obtain at a 95% confidence level. 
As = Q.7Q^q\j and a = -0.09^-^^, a parameter range within 
which the gCg could be a candidate for UDME. However, the 
standard Chaplygin gas with a = 1 is ruled out as a UDME by 
the data at a 99% level. The plan of the paper is as follows. In 
the next section, we provide a brief summary of the gCg and ba- 
sic equations relevant to our work. Constraints from the X-ray 
gas mass fractions in galaxy clusters are discussed in section 3. 
In section 4 we discuss the bounds imposed by the dimension- 
less coordinate distances to SNela and FRIIb radio galaxies. 
Finally, we present a combined analysis, our conclusion and 
discussion in section 5. 

2. The generalized Chaplygin gas: basic equations 

We consider a flat universe that contains only baryonic mat- 
ter and the gCg (we ignore the radiation components in the 
universe that are not important for the cosmological tests con- 
sidered in this work). Then the Friedmann equation is sim- 
ply given by = {^nG/3){pb + PgCg)- Both of the baryonic 
matter and the gCg components satisfy the relativistic energy- 
momentum conservation equation, p+3^(p+p) = 0, where a is 
the scale factor of the universe and '•' stands for the derivative 
relative to cosmic time. From p - for the baryonic matter and 
the equation of state of Eq.(l) for the gCg component, we have 

Pb = Pboa"'; PgCg = Pgcgo [As + (1 - A,)fl-'*'^"')^ (2) 

where pbo and PgCgO are the energy densities of the baryonic 
matter and the gCg at present respectively, and A., = A/p^^^ 
is a substitution of the parameter A. The scale factor is re- 
lated to the observable redshift as a = 1/(1 + z). Now we 
evaluate the dimensionless coordinate distance. v(zL the an- 



Z)^(z), as functions of redshift z as well as the parameters of 
the model. The three distances are simply related to each other 
by = (1 + zfD'^ = (c/Ho)(l + z)y(z). We define the red- 
shift dependence of the Hubble paramter H as H{z) = HqE{z), 
where Hq - lOO/zkms^'Mpc"' is the present Hubble constant. 
The HST key project result is = 0.72 + 0.08 (Freedman et al. 
2001). Parametrizing the model as (A„ or), we get E{z) function 
as (Bento et al.2003a,b; Cunha et al. 2004; Alcaniz et al. 2003) 

E\z\As, a) = Ob(l + z? + (1 - Qb) • 

(a, + (1-A,)(1+z)3(i^«>)^ (3) 

where Q.\, is the density parameter of the baryonic matter com- 
ponent. The observed abundances of light elements together 
with primordial nucleosynthesis give Q-bh^ - 0.0205 + 0.0018 
(O'Meara et al. 2001). Then, it is straightforward to show that 
the distances are given by 

D\z\Ho,As,a) = il+zf-D^iz;Ho,As,a) 
= ^il+z)-yiz;As,a) 

no 

= W^^+^y^w^^^ (4) 
Ho Jo E(z';As,a) 

3. Constraints from the X-ray gas mass fraction of 
galaxy clusters 

As the largest viriaUzed systems in the universe, clusters of 

galaxies provide a fair sample of the matter content of the 
whole universe (White et. al. 1993). A comparison of the gas 
mass fraction of galaxy clusters, /gas = Mgas/Mtot, inferred 
from X-ray observations, with Ob determined by nucleosynthe- 
sis can be used to constrain the density parameter of the uni- 
verse n„ directly (White & Frenk 1991; Fabian 1991; White 
et. al. 1993; White & Fabian 1995; Evrard 1997; Fukugita, 
Hogan & Peebles 1998; Ettori & Fabian 1999). Sasaki (1996) 
and Pen (1997) showed that the /gas data of clusters of galax- 
ies at diff^erent redshifts can also, in principle, be used to con- 
strain other cosmological parameters decribing the geometry 
of the universe. This is based on the fact that the measured 
/gas values for each cluster of galaxies depend on the assumed 
angular diameter distances to the sources as /gas [D'^]^^^. 
The ture, underlying cosmology should be the one which make 
these measured /gas values to be invariant with redshift (Sasaki 
1996; Pen 1997; Allen at al. 2003). However, various uncer- 
tainties in previous measurements have seriously compUcated 
the application of such methods. 

Recently, Allen et al. (2002; 2003) reported precise mea- 
surements of the /gas profiles for 10 relaxed clusters determined 
from the Chandra observational data. Except for Abell 963, the 
/gas profiles of the other 9 clusters appear to have converged or 
be close to converging with a canonical radius f"25oo, which is 
defined as the radius within which the mean mass density is 
2500 times the critical density of the universe at the redshift of 
the cluster (Allen et al. 2002, 2003). The gas mass fraction val- 
ues of these 9 clusters are shown in Figure 1 . With the reduced 
svstematic uncertainties. Allen et al. ^2002: 2003~) successfuUv 
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• Allen et al. (3002) 
□ Allen et al. (2003) 
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Fig. 1. The apparent redshift dependence of the /gas measured 
at r25oo for 9 clusters of galaxies with convergent /ga, pro- 
files. The error bars are the symmetric root-mean- square Icr er- 
rors. The solid circles mark the six clusters studied by Allen et 
al.(2002), while the empty squares mark the other three clusters 
published by Allen et al. (2003). The solid curve corresponds 
our best fit to the gCg model with A, = 0.74, and a = 0.00. 



and Pen (1997) to the data and obtained a tight constraint on 
Q.m and an interesting constraint on cosmological constant. We 
will use this database to constrain the gCg model as a UDME. 
Following Allen et al. (2002), we have the model function as 



(i-H0.i9/ii/2)nf 



^SCDMfc) 



0.5Z^^g(z,.;A„a) 



3/2 



(5) 



where the bias factor b = 0.93 ± 0.05 (Bialek et al. 2001; 
Allen et al. 2003) is a parameter motivated by gas dynam- 
ical simulations, which suggest that the baryon fraction in 
clusters is slightly depressed with respect to the Universe as 
a whole (Bialek et al. 2001). The term represents 
the change in the Hubble parameter from the defaut value of 
Ho = 50kms"^Mpc"' and the ratio Dg^Pj^(z,)/Z)^Pg(z,; Ai,a;) 
accounts for the deviations of the gCg model from the default 
standard cold dark matter (SCDM) cosmology. Note that £1^ 
is the efl'ective matter density parameter (Cunha et al. 2004; 
Makler et al. 2003b), i.e., the coefficient of the term scal- 
ing as (1 -I- z)^ in equation (3) when the gCg behaves like 
dust or equivalently a 1. It is easy to show that, Qj^f = 
Qb + (1 - Qb)(l - A,)'/<'+"'. We should keep in mind that the 
bias factor value for a gCg model might be different from the 
value given above, which leads to a systematic error in this kind 
of analysis. Because b linearly scales the X-ray mass fraction, 
f„=o. in Ea.f5~). lowerine Traisine) it bv ~ 10% would cause the 



Fig. 2. Confidence region plot of the best fit to the /gas of 9 
clusters published by Allen et al. (2002,2003) - see the text 
for a detailed description of the method. The 68% and 95% 
confidence levels in the Aj-a plane are shown in lower shaded 
and lower -i- darker shaded areas respectively. 

A minimization method is used to determine the gCg 
model parameters A ^ and a as follows (Allen et al. 2003) 

l2 



[/gas"^fc'^s,<3^) /gas.oi] 



-I- 



D.bh^ - 0.0205 



0.0018 



h - 0.72 



0.08 



b - 0.93 



0.05 



6) 



where f^°'^iZi',As,a) refers to equation (5), /gas.o; is the mea- 
sured /gas with the defaut SCDM cosmology, and CTf^, is the 
symmetric root-mean-square errors (/ refers to the rth data 
point, with totally 9 data). The summation is over all of the 
observational data points. 

The results of our analysis for the gCg model are displayed 
in Figure 2. We show 68% and 95% confidence level con- 
tours in the (As,a) plane using the lower shaded and the lower 
plus darker shaded areas respectively. The best fit happens at 
As - 0.74 and a - 0.00. Although the data constrain efficiently 
the parameter plane into a narrow strip, the two parameters. A, 
and a, are highly degenerate. This degeneracy can also be seen 
clearly from the relation, (l-Aj)i/(i+'') = (Q^ff-Qb)/(l-i2b). It 
has been shown that the X-ray gas mass fraction is mostly sen- 
sitive to Q„, no matter what the cosmological model is (Allen 
et al. 2002, 2003; Zhu et al. 2004a,b). In our case, a precise 
determination of Q!^ is expected, hence forming a narrow strip 
in the (A a) plane composed of a bundle of curves given by 
(1 - A,)'^*'^"' =const.. In order to determine A, and a respec- 
tively, an independent measurement of A , or a is needed. We 
will show that, in the next section, the dimensionless coordinate 
distances to SNela and FRlIb radio ealaxies are well aoDroori- 
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4. Constraints from the dimensionless coordinate 
distance data 

Motivated by deriving the expansion rate E(z) and the acceler- 
ation rate q{z) of the universe as functions of redshift, Daly and 
Djorgovski (2003) compiled a large database of the dimension- 
less coordinate distance measurements estimated from the ob- 
servations of SNela and FRIIb radio galaxies, and successfully 
applied it for their purpose. We will show this sample provides 
a precise determination of As, and well breaks the degeneracy 
presented in the X-ray gas mass fraction test. 

The database consists in the 54 SNela in the "primary 
fit C" used by Perlmutter et al. (1999), the 37 SNela pub- 
hshed by Riess et al. (1998), the so far highest redshift su- 
pernova 1997ff presented by Reiss et al. (2001), and the 20 
FRIIb radio galaxies studied by Daly and Guerra (2002). The 
authors used the B-band magnitude-redshift relation, iub = 
Mb + 5 log[c(l +z) -yiz)], to determine yiz) for each supernova, 
where Mb = Mb - 5 log Hq + 25 is the "Hubble-constant-free" 
fi-band absolute magnitude at maximum of a SNIa. For the 14 
supernovae that are present in both the Perlmutter et al. (1999) 
and Riess et al. (1998) samples, we will use their average val- 
ues of y with appropriate error bars (see Table 4 of Daly and 
Djorgovski 2003). Therefore we totally have 78 SNela data 
points which are shown as solid circles in Figures 3. The di- 
mensionless coordinate distances of FRIIb radio galaxies were 
estimated through the method proposed by Daly (1994) (see 
also Guerra, Daly, and Wan 2000; Podariu et al. 2003; Daly and 
Djorgovski 2003). We use their values of y for 20 FRIIb radio 
galaxies obtained using the best fit to both the radio galaxy and 
supernova data (see Table 1 of Daly and Djorgovski 2003), that 
are shown as empty squares in Figure 3. 

We determine the model parameters Aj and a by minimiz- 
ing xHA,; a) = Elf 1 bizr, a,, a) - >'o,]^ /o-^, where yizr. A,, a) 
refers to the theoretical prediction from equation (4), joi is 
the observed dimensionless coordinate distances of SNela and 
FRIIb radio galaxies, and cr, is the uncertainty. 

Figure 4 displays the results of our analysis for the gCg 
model. We show 68% and 95% confidence level contours in the 
(A 5, a) plane using the lower shaded and the lower plus darker 
shaded areas respectively. The best fit happens at Aj = 0.74 
and a = 0.32. It is clear from the figure, that the dimensionless 
coordinate distance test alone constrains A , well into a narrow 
range, but limits a weakly. However, it is just appropriate for 
our purpose, to break the degeneracy presented in the X-ray 
gas mass fraction test of last section. As we shall see in Sec. 5, 
when we combine these two tests, we could get very stringent 
constraints on both Aj and a, hence test the gCg as a UDME 
scenario efficiently. 

5. Combined analysis, conclusion and discussion 

Figure 5 displays the results of our combined analysis of the 
constraints from the X-ray gas mass fractions of galaxy clusters 
and the dimensionless coordinate distances to SNela and FRIIb 
radio galaxies. We show 68%, 95% and 99% confidence level 
contours in the (A,, a) olane. The best fit hanoens at Ac = 0.70 



A^=0.74, a=0.3l3 

Aj.= 1.00 (A dominated) 

Aj|=0.00 (Einstein-de Sitter) 
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Fig. 3. Dimensionless coordinate distances y(z) as a function 
of logz for 78 type la supernovae and 20 FRIIb radio galaxies. 
The sohd circles mark the SNela, while the empty squares mark 
the FRIIb radio galaxies. The solid curve corresponds to our 
best fit to the total 98 data points with A , = 0.74, a - 0.32. The 
database are taken from Daly and Djorgovski (2003). 



eracy between A , and a in the X-ray mass fraction test, and the 
dimensionless coordinate distance test is sensitive to A, only, 
a combination of the two data sets gives at a 95% confidence 
level that A,, = 0.70+5j [^ and a = -0.09+°-^;^, a very stringent 
constraint on the gCg. These are the parameter ranges of the 
gCg permitted by the data as a candidate of UDME, which 
is consistent within the errors with the standard dark matter 
+ dark energy scenario, i.e., the case of a = 0. Particularly, 
the standard Chaplygin gas with a = 1 is ruled out as a fea- 
sible UDME by the data at a 99% confidence level. Using the 
CMBR power spectrum measurements from BOOMERANG 
(de Bernardis et al. 2002) and Archeops (Benoit et al. 2003), 
together with the SNela constraints, Bento et al. (2003a) ob- 
tained, 0.74 < As < 0.85, and a < 0.6, which is comparable 
with our results. 

More recently, Bertolami et al. (2004) analyzed the gCg 
model in the light of the latest SNela data (Tonry et al. 2003, 
Barris et al. 2004). They considered both the flat and non- 
flat models. For the flat case, their best fit values for [A,, a] 
are given by [0.79, 0.999] and [0.936, 3.75] wifli and with- 
out the constraint a > I respectively. Particularly, up to 
68% confidence level, the a = 0, i.e., the ACDM case, is 
clearly excluded, though it is consistent at 95% confidence 
level (Bertolami et al. 2004). The authors considered the sce- 
nario in which the gCg unified all matter and energy compo- 
nents, while in our analysis, only dark matter and dark energy 
are unified as the gCg. This might be one factor responsible 
for the difl'erence between their results and ours. Another even 
more imoortant factor is we make heavv use of the X-rav eas 
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Fig. 4. Confidence region plot of the best fit to the dimension- 
less coordinate distances to 78 SNela and 20 FRIIb radio galax- 
ies compiled by Daly and Djorgovski (2003). The 68% and 
95% confidence levels in the A,-a plane are shown in lower 
shaded and lower -i- darker shaded areas respectively. 



analysis depends on the assumption that the /gas values should 
be invariant with redshift, which has been criticised by a minor- 
ity of workers in the field. For example, a recent comparison of 
distant clusters observed by XMM-Newton and Chandra satel- 
lites with available local cluster samples indicate a possible 
evolution of the M-T relation with redshift, i.e., the standard 
paradigm on cluster gas physics need to be revised (Vauclair et 
al. 2003). We should keep this point in mind when we compare 
the results mentioned above. 

Besides various dark energy models (see, e.g., Peeble and 
Ratra 2003), several possible mechanisms without any DE 
component have been also proposed for acceleration of the uni- 
verse, such as brane world cosmologies (Randall and Sundrum 
1999a,b; Alcaniz, Jain and Dev 2002; Deffayet, Dvali and 
Gabadadze 2002; Jain, Dev and Alcaniz 2002, 2003), and 
Cardassian expansion model (Freese and Lewis 2002; Zhu and 
Fujimoto 2002,2003,2004). However, it must be more interest- 
ing if a UDME can be found in which a single fluid plays the 
role of both CDM and DE. The generalized Chaplygin gas is 
such a intriguing candidate, which deserves to explore its var- 
ious observational effects (Kamenshchik et al. 2001; Bento et 
al. 2002;2003a,b; Bilic et al 2002; Bean and Dore 2003; Cunha, 
Alcaniz and Lima 2004; Makler et al. 2003b). In this paper we 
have focused our attention on two observables, the X-ray gas 
mass fraction and the dimensionless coordinate distance. We 
have shown that stringent constraints on the parameters A j and 
a, that completely characterize the scenario, can be obtained 
from the combination analysis of the X-ray mass fractions of 
galaxy clusters and the dimensionless coordinate distances to 
SNela and FRIlb radio ealaxies. It is natuallv hooeful that. 



Fig. 5. Confidence region plot of the best fit from a combined 
analysis for the dimensionless coordinate distances to 78 SNela 
and 20 FRIIb radio galaxies (Daly and Djorgovski 2003) and 
the X-ray gas mass fractions of 9 clusters (Allen et al. 2002, 
2003). The 68%, 95% and 99% confidence levels in the A,-a 
plane are shown in white, white + lower shaded and white -I- 
lower and darker shaded areas respectively. 

various cosmological observations, one could show clearly if 
this scenario of UDME constitutes a feasible description of our 
universe. 
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